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Abstract 



The development of a Modified Perturbation Theory for QCD, introduced in previous 
works, is continued. The gluon propagator is modified as consequence of a soft gluon 
pairs condensate in the vacuum. The modified Feynman rules for a — 1 are shown, and 
some physical magnitudes calculated with them. The mean value of G 2 , gluon masses and 
the effective potential are calculated up to the g 2 order, improving previous calculations. 
In connection with the gluon self-energy it follows that the gluonic mass shell becomes 
tachyonic in the considered approximation. The constituent quarks masses, produced 
by the influence of the condensate, are also calculated. Results of the order of 1/3 of 
the nucleon mass, are obtained for the constituent masses of the up and down quarks. 
In addition, the predicted flavour dependence of the calculated quarks masses turns out 
to be the appropriate to reproduce the spectrum of the ground states within the various 
groups of hadronic resonances, through the simple addition of their constituent quark 
masses. Finally, the generalization of the theory for arbitrary values of the a gauge 
parameter is explored. 
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Chapter 1 
Introduction 



The development of Quantum Chromodynamics (QCD) has been one of the greatest 
achievements of the High Energy Theoretical Physics in the latest thirty years. As it is 
well known the smallest of the coupling constant at high momenta (asymptotic freedom), 
made possible the development of the theory through the usual perturbative language. 
That situation strongly simplifies the study of such phenomena and the calculated re- 
sults are in very good agreement with the experimental data. However, this so called 
perturbative QCD (PQCD) is far from being able to furnish even a rough description 
for the relevant physics at low energies. The solution of this situation is currently one of 
the main challenges of the Theoretical Physics. 

A relevant phenomenon related with QCD is the color confinement. For the time 
being, there are strong reasons to believe that the relation between the basic quantities 
in QCD, the colored gluon and quarks fields, and the real world characterized by a whole 
variety of colorless interacting mesons and baryons, can be understood by solving the 
confinement problem. The basic picture is that a theory in which SU(3) is an exact 
symmetry, the fundamental fields (quarks and gluons) cannot be associated to physical 
states and then the true physical states consist only of colorless mesons and baryons. 

The above cited limitations of the PQCD means in particular that the usual Fock 
space vacuum for the non-interacting theory, is unable to predicts even approximately 
the real ground state of the QCD (I|]-[[§]. This is at contrast with the case in QED where 
the standard perturbation theory, based in the empty of the Fock space, gives a more 
than good concordance with the experimental results. Then, the nature of the vacuum 
structure is one of the main problems to be clarified and naturally its solution is closely 
related with the color confinement effect. A review of the various models considered to 
investigate the confinement and the vacuum structure in QCD can be found in Refs. 

!,!• 
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Another important problem for QCD is the one related with quark masses. There 
are no free quarks, so it is necessary to take care with the meaning of quark masses. 
They should be understood as mass parameters appearing in the QCD Lagrangian 0- 
(called current quark masses) and their origins have generated polemic for long time 
HlOU . The theoretical ratios between the current quark masses and their absolute values 
can be found with the standard methods of current algebra. Their experimental val- 
ues are determined indirectly and the most recently reports for them will be used 
in the present work. On the other hand the masses considered in phenomenological 
hadron models are called constituent quark masses, which values are higher than the 
ones of current masses, and it is supposed that such values come from non-perturbative 
effects, like gluonic condensate effects for example. Other properties of quarks, like 
magnetic moments, couplings, and electro-weak form factors are not influenced by this 
non-perturbative effects |fP2| . 

In the former works |Tj|-[[nj an attempt to construct a modified perturbation ex- 
pansion for QCD, able to predict some low energy properties of this theory, have been 
considered. The proposed modified expansion [ 13 1 conserves the color £77(3) and Lorentz 
symmetries, and was implemented in order to solve the symmetry limitations of the ear- 
lier chromomagnetic field models J^j- 1| which inspired the search. The similar lack of 
manifest Lorentz invariance had also a modification of the Feynman rules, attempting 
to include gluon condensation, advanced later in which a delta function term, for 
k < A qc d, was considered in addition to the perturbative piece. The expansion proposed 
in a previous work ||13|| , considers a change of the gluon propagator in a term associ- 
ated to a condensation of zero momentum gluons. This change of the usual rules had 
the interesting property of producing a non vanishing value for the gluon condensation 
parameter (g 2 G 2 ) [L7|] already in the first approximations. In addition, in Ref. [H| a 
non vanishing value for the effective mass of the gluons was obtained at the one loop 
level. Finally, a perturbative evaluation of the effective potential as a function of the 
condensate parameter indicated that the condensation is spontaneously generated from 
a zero condensate state. 



At this point it is important to remark that earlier works |TB, |19[ introduced a pure 
delta function at zero momentum as gluon propagator, searching for a model of meson 
resonances. In these cases the constants introduced, multiplying the delta function, were 



different to the one fixed in our previous work 1T4L O], which lead to also different results. 



In particular, the singularity structure of the quark propagator have no pole on the real 
p 2 axis [B1J at difference with our results. 
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Justifying the applicability of the Feynman expansion introduced in Ref. [|T^] was 
advanced later |L4|, |15| . This was in need because after modifying the propagator, it 
was unknown whether or not the initial state generating the expansion was a physical 



state of the theory. In Ref. [14, 15|, considering the operational formulation of QCD 



developed by Kugo and Ojima p0[ 2lJ| , was possible to find a zero momentum gluon 



condensate state able to generate through the Wick expansion the sort of propagators 
considered in Ref. [fT3|| . It was checked that this condensate satisfies the physical state 
conditions imposed by the BRST formalism. The discussion allowed also to a more 
precise characterization of the class of changes admitted in the diagramatic expansion 
by the physical state conditions on the initial state. Specifically it was found that the 
C parameter, describing the gluon condensation, must be real and positive. It should 
be remarked that an analogous proposal for a zero mode gluon condensate state was 
independently advanced in Ref. |22| , in an attempt to give foundation to the propagator 
proposed by Munczek and Nemirovski. 

In the present work the development of the Modified Perturbation Theory ||il3f-||T5 



is continued, and the objectives proposed are: a) To state clearly the details of the 
considered perturbative expansion for a = 1. b) To realize an improved evaluation of 
quantities calculated in Ref. fli~3 |, that is, the mean value of G 2 , the gluon mass, and the 
effective potential, by also imposing the restriction that appeared on the parameter C 
describing the condensation |14], pL5f] . c) To apply the results of the calculation of (G 2 ) 



to evaluate the effects of the condensate on the quark masses, after fixing the parameter 
g 2 C to reproduce the accepted value for {g 2 G 2 ). d) Finally, to explore the generalization 
of the theory for arbitrary values of the gauge parameter a. 

The evaluations presented here are done in a special approximation. It can be defined 
as the one allowed to be considered with precision by the fact that in the present state 
of the analysis the renormalization of the fields and coupling constants has not being 
introduced yet. Therefore, the calculations of the (G 2 ), the quark and gluon masses 
will be searched not only up to the order g 2 but also with the additional restriction 
of disregarding all the terms not involving the condensate contribution. To take into 
account those terms will need the introduction of the renormalization procedure up 
to two loops, because the one loop divergence are unchanged by the use of the new 
propagator. This is out of the scope of the present exploration. It is worth noting that 
the just defined approximation scheme resembles a sort of quasi-classical limit in which 
the condensate would play the role of a macroscopic quantum field. The correctness 
of such a picture is suggested by the initial motivation of the analysis, as directed to 
construct a covariant version of the chromomagnetic models |13| . 
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The exposition will be organized as follows: Chapter 2 is divided in two sections. In 
the first one a review of the previous work [jM], [15j is done. In the second section the 
Feynman Rules for the considered perturbative expansion (a = 1) are stated. Chapter 
3 is divided in four sections. In the first one the mean value of G 2 is calculated and the 
free parameter g 2 C fixed from the accepted value of (g 2 G 2 ). In the second and third 
sections the gluon and quark masses respectively are determined. In the four section the 
effective potential is evaluate as a function of (G 2 ). Chapter 4 is devoted to explore the 
generalization of the theory for arbitrary values of the gauge parameter a. Finally, two 
appendixes are introduced. In the first one it is considered the exposition of the diagrams 
calculated for determining the mean value of G 2 . The second one was introduced for an 
analysis of the most elaborated parts in the calculation of longitudinal and scalar modes 
contribution to the gluon propagator modification, considered in Chapter 4. 
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Chapter 2 

Modified Feynman Rules 



The Modified Feynman rules, advanced in previous works, are stated. The gluon prop- 
agator is modified as consequence of a soft gluon pairs condensate in the QCD vacuum, 
this fact is represented by a delta function in momentum representation. In the first 
section the previous work is shortly reviewed and in the second section the Feynman 
rules exposed. 

2.1 Review of "About an Alternative Vacuum State 
for Perturbative QCD" 

In the previous work [TJ, 13] the construction of a modified vacuum state for perturbative 
QCD was considered. By using the operational formulation for the Quantum Gauge Field 
Theory developed by Kugo and Ojima, a particular state vector for QCD in the non- 
interacting limit, that obeys the BRST physical state condition, was proposed. The 
general motivation for searching this wave function was to give a foundation to the 



perturbative expansion proposed in Ref. O]. The proposed state had the form 



cxp 



ACL+ ACL+ , 

^0,1 ^O,! + 



- 48+ Aa+ 



+ B a 



- AL,a+ _| • _o+ an 
S± T ICq Cq 



0). 



(2.1) 



The state (|2.1| ) is a coherent superposition of zero momentum gluon pairs, expressed 



i- 



Aa+ A L ' a + 
^0,2 5 ^0 ' 



ghosts Cq, C\ 



in terms of the transverse and longitudinal gluon Aq 
and Nakanishi-Lautrup Bq + creation operators of the free theory |20| , 
is colorless because the contracted color index a. Eq. (|2.1|) should be understood in a 



This state 
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sense of limit. The operators are supposed to be evaluated at the smallest non- vanishing 
null momentum p resulting from quantizing the QCD free fields with periodic boundary 
conditions in a large box of volume V — *■ oo. The coefficients of the products of creation 
operators are slightly different from their shown numerical values before the infinity 
volume limit is taken. They reach to the values in expression (|2~I1) when V -> oo (p -> 0) 
in the way determined previously flTJL 



For the construction of the state (|2.1]j , the BRST physical state conditions 

Qb\$) = 0, 

Q c | $) = 0, (2.2) 

were required, in which Q B and Qc are the charges associated to the BRST symmetry and 
ghost number conservation. The operational quantization was considered in a scheme of 



the Gupta-Bleuler type, in which all the gluons are treated on the same footing |20|, 
This corresponds with the selection of the gauge parameter a = 1 in the standard 
functional approach. In justifying proceeding in such a way, the central assumption was 
that under the adiabatic connection of the color interaction, the evolution won't bring 
the state out of the physical subspace at any stage of the connection. Thus, the final 
state will be also a physical one for the interacting theory and the associated Feynman 
expansion should have a physical meaning. It is clear that this approach left out the 
question of the construction of a gauge parameter independent formulation of the theory. 
However, as a minimal logical ground for the physical relevance of the predictions was 
given, the generalization of the analysis for and arbitrary a was postponed. In Chapter 
4 of the present work this generalization is explored. 

The gluon propagator corresponding to the modified vacuum state (|2.1|) was calcu- 



lated [15[ through the expression for the generating functional of the free particle 



Green functions Z (J) pjj] . 

($ | exp I % I d 4 x J2 JM,a ( x ) A 7 ( x ) \ ex P y f d * x Yl JAt,a ( x ) A T (•'' > ) I ^ 

I a=l,..,8 J L ^ o=l,. .,8 

xexpT- fd 4 xd i yJ^(x)D; b u (x-y)r' b (y) 

I a,6=l,..8^ 

= ex p u E / d4xdi y J "' a ( x ) v a > - y) JU ' b (y) } . ( 2 - 3 ) 



a,6=l,..8 ' 



where the usual gluon propagator is denoted by D ab v {x — y) and D ab v {x — y) is the modified 
one. 
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The result obtained for the modified propagator in the momentum representation 

was, 



D± (k) 



in/ 



T2 -iC5(k) 



(2.4) 



( |2.4j ), C is the parameter associated with the gluon condensate, and was deter- 

m 



In Eq 

mined to be real and non-negative 

The modification to the standard free ghost propagator, introduced by the proposed 
initial state, was also calculated. Moreover, after considering the free parameter in the 
proposed trial state as real fl2.1|), the ghost propagator was not modified. 



2.2 The Generating Functional and the Modified 
Feynman Rules, for a = 1 

In the present section the modified Feynman rules for the theory are exposed. The 
generating functional of the Green functions, can be written as 



/ D (A, c, c, *) exp {i f d A x (C + J^A^ + c a rf + ffc a + + f } 



J D (A, c, c, exp {i J d 4 xC} 



(2.5) 



where sources have been introduced for all the fields in the usual manner, and the effective 



Lagrangian is given by p0| , |2T 

C 



Gh 



Cc 



C 



Gh 



-id^ (x) Df (x) c b (x) 



C Q = ¥ (x) (i^D^ - m Q S ij ) & 



x) 



(2.6) 

(2.7) 

(2.8) 
(2.9) 



The sum over the six quark flavors will be omitted everywhere in order to simplify 
the exposition, no confusion should arise for it. The gluon field intensity has the usual 
form 



F; v (x) = d,Al (x) - d v Al (x) + gr hc A\ (x) A° (x) 
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D^, D l i are the covariant derivatives in the adjoint and fundamental representations 
respectively, for the SU (3) group 



Df (x) 



d,5 ab - gf abc A^ (x) 
d^5 ij - igT ij < a A«. 



It should be recalled that as argued in |14|, [15] the physical predictions for the value 
a = 1 of the gauge parameter should have physical meaning whenever the adiabatic 
connection of the interaction does not lead the evolving state out the physical space. 
Therefore here this value of the parameter will be fixed for this Chapter and the next 
one. The generalization is advanced in Chapter 4. 

Then in Q2.5p , the standard procedure of extracting out of the functional integral 
the exponential of terms higher than second order in the fields (vertexes) can be used. 
That is, thanks to the recourse of exchanging the fields by functional derivatives of 
the corresponding sources, and also calculating the remaining Gaussian integral, the 



following expression for the generating functional is obtained [24 



Z [J,V>V>€>t] = ^exp 



$abc & 



5 5 



+ 



m 3 5j a sj b s j c 
s s 



+ 



oG 
D abcd 



5 5 5 5 



4!i 4 5J a 5J b 5J c 5J d 
5 5 



S ?aj 5 



2\i 3 5rj r 5J a 5 (—rj 8 ) 



5Zi 5J a 5 (-£ 




with 



Zo [J, V, V, £, £] = exp 



(2.11) 



and the normalization factor is 

6 



N~ 



exp 



+ 



oG 
°abc 



5 5 



Sabcd $ 



5 5 5 



oGh r 
'-'ras u 



3\i 3 5J a 5J b 5J C 
5 5 



+ 



4!i 4 5J a 5J b 5J c 5J d 
S?aj 5 5 5 



2li 3 5fj r 6 J a 5 {-7] S ) i 3 5& SJa 5 (-&• 



Zo [J,V> V> 



J,f},v,U=o 



where use have been made of the DeWitt compact notation [25], in which a Latin letter (a, b, ...), 
for a field symbolizes its space-time coordinates as well as all its internal quantum numbers. 
The same index in a source or a tensor symbolizes the same set of variables of the kind of fields 
associated with this specific index. For example A a = A^ a (x a ), c a = c a (x a ) and Vl/ 1 = ^(xi). 
Such a procedure was useful for the calculations and should not create confusion. As usual, 
repeated indexes represent the corresponding space time integrals and contracted, Lorentz, 
spinor or color components. 



9 



The following definitions have been used in ( [2. 10 ): 



^ijk... 



4r^r4r- sa w) fOT « 



G, G7», Q and $ 
for a = G, Gh, Q, 



A, c, c, * 



where G, Gh and G; mean the gluon, ghost and quark parts of the action respectively. 

There is only one main element determining the difference of the usual Perturbative QCD 
expansion and the modified one considered in |13]-[15|. It is related with the form of the gluon 



propagator. As proposed in [13] there is an additional term in the gluon propagator which is 
absent in the standard expansion. In [14, 15], such a term was shown to be a consequence 
of a Wick expansion based in a state constructed by acting on the usual vacuum with an 
exponential of zero momentum pairs of gluon and ghost creation operators. By the side, the 



ghost propagator 14, 15] could remain unmodified if the parameter, free after fixing the form of 



the wave- function, is taken as a real and positive one [14, 15]. Such a selection will be employed 
here. In Sect. 3.2 possible implications of the ghost propagator modification will be discussed. 
The quark propagator is not affected in any way by the gluon condensation as introduced in 



From a more technical point of view the operator quantization employed in |14j, |15j takes the 
ghost fields as satisfying the conjugation properties defined by the Kugo and Ojima approach 
pfj| , 21]. However, at the level of the Feynman diagram expansion the difference with the 



standard procedure is only a change of variables. 

In accordance with above remarks and the results of [13|-|15] the diagram technique defining 
the modified expansion has the following basic elements 



Propagators 



/-its 
^Gh 



G) 



k 2 + 



IS 



iC6(k) 



S'-'i 



k 2 -Me' 

m + P^ln 
(to 2 — p 2 — ie) 
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aG 
D, 



abc 



J, 



abed 



s: 



Gh 



He 



vertexes 

(2vr) 4 5 (k a + k b + kc) (-i) gf abc g^ c (k a - k c )^ + g^ b (k b - k a ) 

+W (k c ~k b ) 

(2vr) 4 5 (k a + k b + k c + k d ) {-g 2 ) (g^ c g^ d - Wc <WJ + 

a 

+ / / i.9 HdHc9 HaUb ~ 9fJ,bfJ-c9fJ,afJ.d) 

(2vr) 4 5 (Av + A; a + fc.) (fc r - fc s ) Ma , 
(27r) 4 <Kk + ** + *i)<#tf7/i.- 



The propagators and vertexes will be represented by, 



.ab 



VJfJ 



Gh 



'abc 




abed 




iaj 



Figure 2.1: Propagator and vertex representations. 

Being defined the Feynman rules under consideration, in the next section they will be 
applied to the calculation of various quantities of physical interest. 
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Chapter 3 



Calculations within the Modified 
Perturbation Theory 



Some magnitudes of physical interest are calculated within the modified perturbation theory, up 
to the g 2 order, and in the approximation described in the Introduction. That is disregarding 
all the terms not including the gluon condensate. The mean value of G 2 (Section 3.1) and the 
gluon mass (Section 3.2) are determined improving previous calculations p3| . In Section 3.3 
the constituent quarks masses are calculated. Finally the effective potential [13] is obtained in 
terms of the mean value of G 2 (Section 3.4). 



3.1 Mean value of G 2 

The calculation of the mean value of G 2 is considered in this section. It completes the former 



evaluation done in [ 13 1 in the tree approximation. The result will be employed afterwards for 
fixing the value of the condensation parameter g 2 C in order to obtain the currently accepted 
value for (g 2 G 2 ). As it was stated in the introduction, the quantities will be exactly calculated 
disregarding all the terms not containing the influence of the condensate. It can be remarked 
that all the terms up to the g 2 order considered here, in which the condensate parameter enters 
and using dimensional regularization, do not need for renormalization. This result is another 
resemblance with a quasi-classical approximation. 

For the calculation of (G 2 ) the following expression will be used 



<0| S g |0> 



/ D (A, c, c, *) S g [A] exp {i J d 4 x (C + J^MJJ + c<V + fj a c a + + } 



/ D (A, c, c, *) exp {i J d A xC] 



(3-1) 
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where S g [A] represents the gluon part of the action in absence of the gauge breaking term 
depending on the gauge parameter a, that is 



S g [A] = -\J d 4 xF- u (x) F^ (x) . (3.2) 
Through the use again of the trick of changing the fields inside the functional integral by 



derivatives in corresponding sources, it follows for (3.1) 



<0| S g |0> 



S 9 ab 8 6 . £&. 5 5 5. S±, 5 5 5 5 



+ 



abc 



+ 



abed 



2i 2 5J a 5J b 3H 3 5J a 5J h 5J C 4!i 4 5J a 5J b 5J C 5J ( 



z [J, V, V, C, f] 



(3.3) 



where Z \j, 77, 77, £, ^1 is given by the expression ( 2.10| ). In the momentum representation S 9 ab 
has the form 



S 9 ab = (2vr) 4 5 (k a + k b ) 5 ah [- [g^kl - ka^ka^)] , 
and it will be represented in the calculations by 

Figure 3.1: Representation of Sab- 
It should be noticed that the contribution of diagrams including any kind of particles (which 
involve the condensate effects), to the vacuum expectation value of S g , needs for the evaluation 
of loop integrals which only can be properly accounted after renormalization. However, in the 
present situation for these contributions to (0| S g |0), it happens that the loop integrals which 
appears are evaluated at zero momentum. Then, the results vanish using the dimensional 
regularization and thus do not depend on the renormalization conditions. 

After computing all the diagrams appearing in the (0| S g |0) expression, with the remarks 
made above, was obtained that the only diagram contributing is the one appearing in the 



Figure [3^. That is, all the diagrams not considered in earlier works have a null contribution. 
The details of the calculation are exposed in the Appendix A. 

<ols g lo>=-loO 

Figure 3.2: Diagram contributing to (0| S g |0). 
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The diagram in Figure 



3.2 



was calculated before (T^, [14|], the result is 



(0\S 9 \0) = -\jd^ 2 -^^, (3.4) 
which, in turns implies for the mean value of G 2 the result 

(g2> - W- (3 - 5> 

Finally, making use of a currently estimated value for (g 2 G 2 ): 

(g 2 G 2 )^ 0.5 (GeV/c 2 ) 4 , 
the parameter g 2 C is determined to have the value: 

g 2 C = 64.94 (GeV/c 2 ) 2 . (3.6) 

The fixing of the condensation parameter allows to investigate the physical predictions 
of the modified expansion in the above described approximation which does not require to 
introduce the renormalization. 

3.2 Gluon Mass 

In the present section it is exposed the calculation for the gluon mass. 

The expression for the effective action of the gluon sector, in the one loop approximation 



or up to g order |j25 |, is 

1 



r G [4> c ] = S G [4> c ] + i: In det G G [0 C ] , (3.7) 



in Feynman diagrams this expression can be written in the form, 



r [#] = • + ! O 

2i 

Figure 3.3: The effective action. 
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The second functional derivative of (3.7), in the gluon fields, defines the two point proper 
Green function or the inverse of the gluon propagator. It can be written in the form 



r£ [<Pc] = S% [<(> c ] + I (S% qj [<p c ] G" + S% q [<f> c ] G%G^S% [<p c ]) . (3.8) 

where the quark and ghost fields do not contribute in the approximation considered. That is, 
their diagrams do not involve the gluon condensate. 



As it was stated in a previous work [13], the interest is centered in the case of zero mean 
values of the gluonic field (as it is required by the Lorentz invariance), so the expression ( p.8| ) 
takes the form, 

[0] = &ij + ^ (Sf pq jG P Q + Sf pq G q G Gq S sr j} , (3-9) 

or diagrammatically 



Figure 3.4: Two point proper Green function. 

Writing the expressions in the momentum representation, the following results are obtained 
for the gluon condensate contributions 

S G = -5^g^ N p\ (3.10) 

aG r vq _ a ^9 2 GN 
(2vr) 

S? pq G%G^Sf rj = i^S^ (-10^ + 4^) . (3.12) 
The sum of the three contributions is 



r%[0;p]=8 a ^ 



2 | WCN \ PjhPvA q PihP^ 

+ (27r) 4 ) [9** p2 ) P p 2 



(3.13) 
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For calculating the previous expressions the following relation, between the structure con- 
stants for a gauge group SU (N) , was employed 



Considering the specific group SU (3), the expression ( 3.13 ) takes the form 



(3.14) 



T%{0;p]=S^ 



2 6^C\ / PjhPvA q PjmPm 



p 



(3.15) 



The eigenvalues of the previous expression define the dispersion relations, and have the 
form 



p 2 = 0, p 2 



where 



m 



G 



(3.16) 



m 2 G 



6g 2 C 



-0.25 (GeV/c 2 ) 2 . 



(3.17) 



Therefore, as the parameter C was defined before as real and positive, in the framework of 



the initial state considered [14, |i5[ , it follows that the transverse gluon mass correction becomes 
tachyonic. The ability of a tachyonic mass in producing improvements in models for the inter 
quark potential has been recently argued in the literature The effect is related with 

the introduction of a linearly rising term as a first correction to the Coulomb potential in the 



massless case 27]. Therefore, the tachyonic result arising here for the gluon mass, appears 
to be of interest for an attempt of deriving from the present framework the already existing 
successful phenomenological bound state models for mesons [28|-|3C]. 

Finally it should be mentioned that if a modification of the ghost propagator is considered, 
then appears a longitudinal mass term for the propagator, as consequence of the diagram in 



Fig. 3J5 and the polarization operator (gluon self-energy) becomes non-transverse. This result 
means a violation of a Ward identity representing the gauge invariance, by the presence of the 
ghost condensate. Thus, it follows from the present results that the selection of the parameter 
for the initial state in [|l4|, 15] was the appropriate one for the fulfillment of a gauge invariance 
condition at this stage. The physical meaning of alternative procedures for constructing the 
initial state need however a clarification. 
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Figure 3.5: Ghost contribution to the gluon self-energy. 



3.3 Quark Masses 



The effects produced by the condensate, on the quark masses, are considered in this section. 
The full inverse propagator for quarks, in terms of the corresponding mass operator, can be 
expressed as [31] 



G Q 1 (p) = (m Q -p^ ll -Z(p)), (3.18) 
where the color index identity matrix has been omitted and E (p) is the mass operator. Its 



expression is represented by the diagram in Figure 3.6. 



il_ i} ^ 



Figure 3.6: Diagram for E (p). 
The exact form of E (p) up to the g 2 order, considering the modified Feynman rules, is 

S (P) = 9 2 C F [ d*k ^{rn Q + (p-kT la)l G%(k)^ 



in which the algebra for the color indexes has been calculated 31], leading to the following 
expression for the constant Cf (SU (N)) 

_N 2 -1 

for the case of interest SU (3), it reduces to Cp = 5. 
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Taking into account the modified gluon propagator and the standard relations for the 7 
matrices (which here are adopted the same as in Ref. |31|) the full one loop mass operator 
expression can be simplified to 

" W = & r &■-*>>> W> ■ (3-20) 



(2n)U (m 2 Q -(p-k) 2 ) V^ 2 



Here, as before, the term not involving the condensate won't be considered. In this case 
the expression for £ (p) reduces to 



g 2 C F C2(2m Q -p a % 
(2vr) 4 ^ m ^_p2^ 

Then the expression for the inverse quark propagator reduces to 



(27r) 4 ( m i - p 2\ 



where it is defined 



Gq 1 (p) = m Q \l + 2 ( M - 1 + M I , I :*>.22 , 

m Q ~ P ) I \ ( m Q ~ P 



m 2 = ^CfC =qa 2 2 

(2vr) 4 V ' ' 



in which the numerical value has been obtained considering (J^ 

The zeros of the determinant associated to the inverse propagator ( |3.22| ) allow to determine 
the effects of the condensate (reflected by the parameter M) , on the effective mass or mass shell 
of quarks. Now is useful to introduce an index s,s = 1...6 for each kind of quark characterized 
by its particular current mass mg s . 

The mass shell obtained in the present case is 

p 2 - ml . = 0, (3.23) 

for m 2 a . are obtained three different analytical expressions arising from the solutions of a cubic 
equation. They are 

1 2 
m 2 = A^-B + m 2 n +-M 2 , 

9s, 1 Vs g ' 

2 1 ' 



—A* + -B + m 2 n + -M 2 + -i>/3 (a* + b) , 
2 2 Qs 3 2 V J ' 

I, = — A% + -B + m% +-M 2 - -iy/3(A* +b] , 
2 2 Qs 3 2 V / 



ys,3 
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in which A, B are given by the following expressions 



A 
B 



-m 2 n M 4 M 6 + —J96m% M 6 + 177m 4 M 8 

6 Qa 27 18 V Qs Qs 



l2m 2 Qs M w , 



±M 4 



aI 



In order to find the dependence of the effective quark masses m 2 g . as functions of the 
current mass parameters m,Q g , an analysis for the various dispersion relations that appeared, 
was done. It follows that there is only one solution having a squared mass being positive for 
arbitrary values of the current quark masses. The existence of this real and positive solution 
for the squared mass is possible because of the real and positive character of C, as was proved 
previously Ji~4| , 15] and fixed through (g 2 G 2 ). The other solutions for the masses become 
complex for certain values of the current masses m,Q a , and there is no a clear understanding of 
their meanings. 

For the purely real solution the value of the effective quark masses m qa = f {rriQ 3 ) as 



functions of mQ s is shown in the Figure 3.7. The graph is plotted for the region mQ a < 2 GeV/c 2 
which contains the current mass values of the u, d, s and c quarks. 




Figure 3.7: Real solution for the quark mass as a function of the Lagrangian mass (masses 
in GeV/c 2 ). 
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As it can be appreciated in the picture, the calculated effective quark masses for light 
flavors (u, d and s) are clearly predicting the values of the quark masses being in use in the 
constituent quark models of hadrons. That is, the light quarks get a weight of near one third 
of the nucleon mass. In table 3d the mass values obtained, for each quark flavor, are shown. 
The reported lower and upper values for the current masses are denoted by and 
respectively fll|]. The correspondingly calculated constituent masses are indicated by rn^ w 
and rnjj h p ] finally rnj^ ed is the evaluated mean value of the constituent masses. 



Table 3.1: Quark mass values (MeV/c 2 ) in presence of the condensate. 



Quarks 


h)xp 
m Low 


tJxp 

m u P 


Th 
m Low 


Th 
m Up 


Th 
" b Med 


Up (u) 


1.5 


5 


334.1 


335.8 


334.9 


Down (d) 


3 


9 


334.8 


337.8 


336.3 


Strange (s) 


60 


170 


361.3 


415.1 


388.2 


Charmed (c) 


1100 


1400 


1172.6 


1457.9 


1315.2 


Bottom (b) 


4100 


4400 


4120.3 


4418.9 


4269.6 


Top (t) 


168600 


179000 


168600.5 


179000.5 


173800.5 



From the global properties of table |3d] it can be observed that the main effect of the gluon 
condensate seems to be dressing the light quarks with a cloud of gluons having a total mass of 
one third of the nucleon mass. Then, the results point in the direction of the mainly glue nature 
of the constituent quark masses of the u, d and s quarks within many baryon resonance. These 
results support the idea that a modified perturbative expansion like the one being considered, 
in which the effect of the gluon condensate has been incorporated, would be able to predict 
with reasonably good approximation properties of the low energy strong interactions. Similar 
values for the constituent quarks masses have been obtained by a different method in [32]. 



In tables |3.2| and 3.3 are shown the masses for baryon and meson ground states [11 



compared with the values obtained for them by a theoretical estimate obtained as the simple 
addition of the calculated constituent masses for each of their known constituent quarks. 
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Table 3.2: Experimental and theoretical masses (MeV/c 2 ) for baryonic resonances in 
their ground states. 



Baryon Exp.Val. 


in. Val. 


KeLiirr. 


p(uud) 


938.27231 


1UU6.Z 


1 .1 


n(udd) 


939.56563 


1UU / .0 


l.Z 


A(uds) 


1115.683 


iUoy.4 


O.U 


E + (uus) 


1189.37 


iujo. i 


11. u 


S u (uds) 


1192.642 


1059.4 


11.2 


S-(dds) 1197.449 


1060.8 


11.4 


H°(uss) 


1314.9 


1111.4 


15.5 


S~(dss) 


1321.32 


1112.7 


15.8 


Q-(sss) 1642.45 


1164.6 


29.1 


A+(udc) 2284.9 


1986.5 


13.1 


2+ (use) 2465.6 


2038.4 


17.3 


H^(dsc) 2470.3 


2039.7 


17.4 


^(ssc) 2704 


2091.6 


22.6 


A^(udb) 5624 


4940.8 


12.2 



Table 3.3: Experimental and theoretical masses (MeV/c 2 ) for vector messon resonances 
in their ground states. 



Meson 


Exp.Val. 


Th.Val. 


Rel.Err. 


,_, ( uu—dd \ 
P\ V2 ) 


770.0 


671.2 


12.8 


„ ( uu+dd \ 
W \ V2 ) 


781.94 


671.2 


14.2 


<j)(ss) j 


1019.413 


776.4 


23.8 


Jji> (IS) (cc) 


3096.88 


2630.5 


15.1 


Y(1S) (bb) 


9460.37 


8539.1 


9.7 
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As it can be appreciated the results obtained for the vector meson and baryon resonances 
are only estimates. However, their values are in good correspondence with the reported exper- 
imental values. This situation is more curious taking into account: 1) The great errors in the 
quark current masses |yj 2) The standard one loop self energy contribution to the propagator 
was not included. The major differences between the experimental values and the theoretical 
calculations are obtained for resonances containing the strange quark. The improvement of 
these results could be introduced from higher order corrections, because the strange quark is 
located between the almost massless up, down and the heavy quarks, then the one loop correc- 
tion could not be enough to obtain the total influence of the condensate on that quark. It can 
be noticed that the validity of the results obtained implies that the binding energy contribution 
to the baryon rest masses is small. This property then could justify a no relativistic description 
of such resonances. 



3.4 Effective Potential 

In the present section is calculated the effective potential as a function of the mean value of G 2 



The approximation considered is the same advanced in a previous work [13], and the result is 
improved through the relation between the mean value of G 2 and the gluon mass obtained in 
previous sections. It is shown that the condensate is spontaneously generate from the standard 
vacuum, as in the Savvidy models Q. 

The expression for the effective potential, considering the presence of a mean field, was 
obtained by R. Jackiw and has the form 

V (0 C ) = V (0 C ) + j.f -0^ lndet [<f> e ; k]) + i (exp (i J d^xQ [</> e , <f>]j ^ , (3.24) 

as it was mentioned before, in the present work is analyzed the case for null mean values 
of gluonic fields {(p c = 0), as it is required by the Lorentz invariance. The first term in 
the previous expression is (G 2 )/4, the second term represents the quantum corrections to the 
effective potential up to one loop approximation, and the last term represents the higher order 
corrections. 



The approximation considered in the present section [13|, is the one in which all the mass 
insertions, consequence of the gluon condensate, are introduced in the one loop term of the 
expression ( [3~24D . That is, it is substituted the 5"?- [0; k] by the T^ - [0; k] calculated in Section 



3.2, and it is like considering T G ^ [0; k] as a modified tree approximation for the second derivative 
of the action under the influence of the condensate. These corrections appear in an exact 
calculation from the expansion of the third term of ( |3.24 ), that won't be taken into account in 



what follows. In this case the expression ( |3.24j ) takes the form 

?2\ 1 r rl^l 



^«G 2 » = ^ + ^ / " lndet (-rg[0;fc]), (3.25) 
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that considering the relation 

lndet (-rg [0; k}) = Trln (-Tg [0; k]) , (3.26) 

can be written as 

^((^ 2 )) = ^ + ^ / ^Tlln(-rg[0;ib]), (3.27) 

In order to calculate the previous integral, the relation ( |3. 15| ) can be written in the form 

T% [0; k] = -5 a ^ [(k 2 - m 2 G ) P T + k 2 P L ] , (3.28) 
where Pt is the transverse projection operator 

PT = 9, m -^, (3-29) 
and Pl is the longitudinal projection operator 

Pl = (3.30) 
These operators satisfy the following relations: 

(P T ) N = P T: for N > (3.31) 

(P L f = P L , foriV>0 (3.32) 

Pt-Pl = 0, (3.33) 

Tr(5 a ^P T ) = 8 x 3 = 24, (3.34) 

Tr(<5 a ^P L ) = 8x1 = 8. (3.35) 

Then it is possible to write 

Trln (-rg [0; k]) = Trln (5 a ^ [(k 2 - m%) P T + k 2 P L \) 

= 24 In (k 2 - m 2 G ) + 8 In (k 2 ) , (3.36) 

and considering that 

j- In (A; 2 ) = const, (3.37) 

(2tt) i 

the attention will be centered in the term 

f d 4 k . » 2 , f d A k ( m 2 r \ 
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For the calculation is defined: 



d 4 k 
(2^) 4 i 



In 1 



m 



G 
k 2 



(3.39) 



It is not easy calculate the integral in the previous expression, because it should be compute 
in the Minkowski space. The solution of this problem is obtained through a 90° rotation in 
the complex plane ko, with this transformation the Minkowski integration is replaced by an 
Euclidean integration. Such a rotation is called Wick rotation and is equivalent to a change of 
variables 



o. 



with Kq real 



then 



d 4 k = id 4 K, k 2 = -K 2 , with K 2 = K 2 + 1? 2 . 



and (3.39) takes the form 



d 4 K ( m 2 r 
In 1H — 



(2^ 



(3.40) 
(3.41) 

(3.42) 



Differentiating the previous expression in m G it is possible to eliminate the logarithm in 
the integral 



dT m 
dm 2 G 



d K 



1 



K 3 dK 



(2n) A K 2 + m 2 G J K 2 + m% 



and considering for the angular integration the expression [31]: 

<ffi _ 2vr 2 
(2^?- (2^' 



it is possible to rewrite (3.43) in the form 



dT m 
dm 2 G 



2tt 2 
(2nf 



K 3 dK 
K 2 + 



in 



G 



dQ 

(2^r 



(3.43) 



(3.44) 



(3.45) 



The momentum integral in ( 3.45| ) is divergent, so it is needed to use dimensional regular- 
ization to calculate it. The area of a L>-dimensions unitary sphere is 



S D - 



D 

2vrT 



(3.46) 



24 



then the integral ( p. 45 ) in D-dimensions has the form (2 



dT m 



S D - 



(27T) 



D 



K D ~ x dK 
K 2 + 



m 



2-E 

G 



m 



G 



(4vr 



,2-e 



(3.47) 



where e = 4 — D and in general 



-n + 
* (n + 1) 



i 



+ *(n + l) +0(e) 



1 1 

1 + - + ... + - -7, 

2 n 

-7 with 7 = 0.5772 



(3.48) 



All the functions in ( 3. 47] ) are expanded for e — > 0, and cutting the pole in e the expression 
( 3.47p takes the form, 



dT„ 



drriQ V 4tt ) 



(3.49) 



and the result for the integral in uiq is 



Tm= 64^ m ^ ln 



G 



const 



(3.50) 



Finally, the expression (3.24) can be written as 



V UG 2 )) = ^ + 12-^m G ln (3.51) 
vx /; 4 64vr 2 G const v ; 

that considering ( |3.17| ) and ([0]) takes the form 

which has a minimum for G 2 different from zero. This implies that the gluon condensate is 
generated spontaneously from the zero condensate vacuum, as for Savvidy models. 
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Chapter 4 

Generalization of the Theory for 
arbitrary values of the a parameter 



The generalization of the proposed theory, for arbitrary values of the gauge parameter a, is 
explored. In Sec. 4.1, the ansatz for the Fock space state that generates the desired perturbative 
expansion is defined. The proof that the state satisfies the physical state conditions imposed 
for the BRST formalism is also analyzed in this section. In Sec. 4.2 the modifications to 
the usual perturbation theory are calculated, and the introduced parameters fixed in such a 
way of making explicit the Lorentz invariance in the propagator and reproducing a previous 
proposition fi"3| |. 



4.1 Modified Initial State 

In the present section it is considered the physical vacuum state able to modify the usual 
perturbative theory in the way proposed previously |l3|], for arbitrary values of the a parameter. 
Some needed results of the Kugo and Ojima quantization, for arbitrary values of the gauge 
parameter, are also pointed out. 

The equations of motion for arbitrary values of a, deduced from the Kugo and Ojima 



Lagrangian [20, 21 1, are 



nA«(x)-(l-a)d,B a (x) = 0, 
d»A a ^(x) + aB a {x) = 0, 
UB a (x) = Uc a (x) = Oc a (x) = 0. (4.1) 

Only the solution obtained for the gluonic field differs from the one obtained for a = 1, 
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considered earlier [TL4L uM , and takes the form Ml u\ 



k V=l>2 

+h.c, 



x 



where the wave packets in (|4.2| ) are taken as 

, , 1 



9k [x) 



exp (—ikx) , kx = koxo — kx, ko 



^/Wk~ Q 

/fc,/i 0»0 = Tyfp^^ exp H fca • 
The polarization vectors (fe), e/, jAi (A), (fe) are defined by 



k ■ ^ (k) = 0, e a (k) = 0, e° (k) ■ ? (k) = 5 aT , for a, r = 1, 2 



—1 



,k 



(4.2) 



(4.3) 



(4.4) 



The commutation relations between the creation and annihilation operators are given by 



At 

k,a 
k 


A a '+ 

k>,a> 
^'^^ 




A L,a' + 
k' 





Bi' + 
k' 




c a + 
h' 





c a + 
k> 





BZ 

k 



















^ aa 'hk 













-iS aa '5 rk , 





The integro-differential operator 


D {1/2) is 


defined by 







(4.5) 



DfJ } 2) = \/2{V 2 )-\x d - 1/2), 



and works as an "inverse" of the d'Alembertian for simple pole functions [34] 



(4.6) 



UD^ f (x) = f {x) if Df(x) = 0. 



(4.7) 



27 



The value of d^D^l 2 ^ gk{x) in (4.2) is determined with the use of the wave packets (4.3) 



and polarization vectors ([4.4[) considered. The result obtained is 



d.D^guix) 



d^/2{V 2 )- l {x d -l/2) 



1 



exp (— ikx) 



A (ko,x ) fkL.u, (a?) + Bf, (k, x) 



where A (ko,Xo) and B^ (k, x) are defined by: 

A(k ,x ) = — 
13^ (k,x) 



(ik x + 1/2) , 



ikoSuo , . 



In similar way for c^-D^Y 9% ( x ) ^ s obtained: 



(4.8) 



d^D^gUx) = A* (k , x ) Hl^ (*) + (fc, x) . 



(4.9) 



Now it is possible to define the proposition for the QCD vacuum state. The considered 
one is a generalization of the previously proposed [|i~4"| , |l5| , in which the constant accompanying 
the non-physical sector has allowed to be gauge dependent in order to obtain a covariant 
gluon propagator. This state is physically equivalent to the one proposed in the previous work 
[14, 15], because its physical part (transverse gluons part) was not modified. The modifications 
were introduced only in the non-physical sector, as usual for this sector, to make explicit the 
covariance. 

The new vacuum state is defined by the expression 



[ a=i Pi,\pz\<p 

In this state the gluon pairs are created with very small momenta pi, but different from zero, 
with the objective of avoiding the difficulties appearing if they are considered with zero mo- 
mentum directly. The sum over pi for \pi\ < P is introduced to avoid any preferential direction 
in the space when the limit P — > is considered, taking into account that the volume V —* oo. 
Is interesting to note that this sum was not introduced in a previous work (|l4| , ftfj, for the 



E 

<r=l,2 



cap) 



A a+ A a+ + c / m f B o+^,<*+ + , 
Pi,a p,,cr -5 V ' > \ p, Pi Pi -j 



a+ „a+ 

Pi 



0). 
(4.10) 
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a = 1 analysis, because the operator D^t( was not involved and the results were independent 
from the introduced auxiliary momentum direction. 

In the analysis following the work will be explicitly done with a generic color and a generic 
momentum pi, for simplifying the exposition. This can be done because in the free theory all 
different colors and momenta can be worked out independently, thanks to the commutation 
relations defined (|4.5|). At the necessary point of the analysis all the contributions will be 
included. 

The proof that the proposed state (4.10) satisfies the required physical state conditions, 



Qb I *) = 0, 

Qc\V) = 0, (4.11) 

is exactly the same advanced in the previous work, because the constant C3 (a, P) does not 
make any change. So this state is a physical state of the theory and it is possible to state that 
the non-physical sector of | ^) is undetectable in the physical world and it can be neglect in 
physical calculations, like in the calculation of the norm for example. 

Then for the calculation of the norm a similar result, to the one in [jl4|, |l5|| , is obtained 

iv=<*i*>= n n n 

CT =l,2o=lp i ,|ft|<P 

The normalized physical vacuum state is then defined by 



P 



1 2m 



(2m)! 



for \C a {P)\<\ 



m=0 



= -= I m. (4.12) 

y/N 



4.2 Modified Perturbation Theory 

In the present section the modification to the usual perturbation theory is calculated. It is 
introduced through the modification of the usual vacuum state (empty of the Fock space), by 
the state proposed in the previous section ( 4.10| ) . The analysis is generalized for arbitrary values 



of the gauge parameter a. At the end a brief analysis of the a generalization consequences, on 
the previous calculations for (G 2 ), quark and gluon masses, is done. 
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As it was pointed out in a previous work ]i"4] , [l5| , the modification to the gluon generating 
functional of the free-particle Green functions (and to the propagator also) is determined by: 

($ | exp i i j d 4 x Jfl,a ( x ) A T ( x ) \ ex P \ 1 [ d * x E J ^ ( x ) A t ( a 

[ ^ a=l,..8 J [ J a=l,..8 



Considering the expressions ( |4.2| 
operators in Q4.13j ) have the form 



tf). (4.13) 

and (|4.9| ), the gluon annihilation and creation 



K + (?) 



E E^L/^)+4' a W*)+ 



, (7=1,2 



+52 [f ktSifl (x) + (1 - a) (A (A*,, x ) fk,L,n (x) + fl„ (k, x))] 

E ( E A TJ^ (*) + 4 0+ ^*^ (*) + • 

icr=l,2 

+ (!-«) (*d,a;o)/ fc *r lM (ar) + S*(*;,x))] ) .(1.1 I) 



As it can be noticed in the expressions ( [4,10| ) and ( 4.14j) the transverse and ghost terms 
are not modified with the generalization of the theory for arbitrary values of a. Then their 



contributions are the same calculated previously [14, |lq| . It is important to remark that the 



calculations are done decomposing the expression (4.13) in products of each mode contribution, 
thanks to commutation relations ( f4.5[ ). 

The result obtained for the transverse mode contributions, for \C a (P)\ < 1, is |l4| , |D| 

(° I ex P ( E ^^ A k« A k^j ex P |< / d'xJ^{x) e« ( ,)} 

x exp i i f d A xJ^ a (x) Y, A kJL, (*) \ ex P ( E 
[ <7=l,2 J \(7=l,: 



— — yl a+ A a+ 

2 Pii " Pii 17 



exp < 



E CCO 5 

(7=1,2 



C a (P) + C*(P) + 2\C a (P)\- 



2 ll-\C a (P) 



(4.15) 



where the normalization factor has been canceled, and it should be reminded that the momenta 
Pi are very small and the sources are located in a space finite region. The following simplified 
notation was also introduced fbil, |l5| 



ja 



d\ 



-.J^ a (x) e a „ ( Pi ) . 
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The calculation of longitudinal and scalar mode contributions is in this occasion more 
elaborated, and a brief exposition of it could be found in the Appendix B, the result obtained 
for \C 3 (a,P)\ < 1 is 



exp < 



d A xd A y 



2Vp 



■iO 



j»> a (x) r> a (y) 



C 3 (a,P) + C* 3 (a,P) + 2\C 3 (a,P 



x es,M (Pi) £l,v (Pi) + 



(!-<*) 

m\ 2 



l-\C 3 (a,P) 



[ez,!* (Pi) + i^PioSfjo] e L)V (p. 




(4.16) 



the normalization factor has been canceled. 

Now substituting the expressions pill ) and pl6| ) in ( glf ). Performing some algebraic 
manipulations, keeping in mind the properties of the polarization vectors defined (4.4), consid- 
ering C\ (P) = C2 (P), and finally introducing the contributions for all momenta pi for \pi\ < P, 
the following result is obtained 



exp 



C 1 (P) + C{{P) + 2\C 1 {P)[ 



+ 



+ 



2(l-|d(P)| 2 

C 3 (a, P) + C* 3 (a, P) + 2\C 3 (a, P)\ 2 C x (P) + C{ (P) + 2 |d (P)| 2 \ />„,/>„, 



l-\C 3 (a,P) 



C 3 (a,P)+CI(a,P) + 2\C 3 (a,P)\ 2 \ (1 - a) 



l-\C 3 (a,P)[ 



1 - |d (p) 



(Pi/iPii/ - 2p i0 5 f _ l op il/ ) 



2\A 



(4.17) 



It is interesting to note at this point that the combinations of C\ (P) and C 3 (a,P), in the ex- 
pression ( |4.17| ), are real. So even if these parameters are complex the modifications introduced 
by them are real. Then, after knowing that complex parameters essentially produce the same 
result as the real ones, they will be considered as real in what follows. 

In the expression Q4.17 ) the addition is changed by an integration, considering that V — ► 00, 
in accordance with the relation 



- Y, 



V 

Pi,\fi\<P 

and the integrals are calculated. 



2tt 



(2ttV 



dp p" 



J sin 6d6 j dip, 



(4.18) 
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After that, the parameters introduced were considered behaving in following form, for 
P ~ 0, 



Ci(P) 
C 3 (a,P) 



~ 1 - ^-P 2 , Ci > 0, 



c 3 («) o2 



P', C 3 (a) >0. 



The result obtained for the expression ( 4.17 ), in the limit P — > 0, is then 



(4.19) 



exp 



<FxctyJ^ a (x) J^'" (y) 



+ 



1 



1 



1 



(2tt) 2 Ci 3(2^) 2 V^3(a) C] 
(1-a 



(2vr) 2 6C 3 (a 



(5^ + 2 <^mo^o) 



(4.20) 



The constant C3 (a) is fixed in order to cancel the terms proportional to SuoSvo, to make 
explicit the Lorentz covariance. That's why this non-physical part was introduced in the 
vacuum state. Then C3 (a) is determined by the expression 

(1 + a) 

L-3 (a) = 7: t^l, 



and ( 4.20 ) takes the form 



exp 



cPxtfyJ** (x) J u ' a (y) 



(2^) 2 Ci 



(4.21) 



In the expression (4.21) is important to note that the term in the exponential is real and non- 
negative Q4.19p . And that it is independent of the gauge parameter, which is a very interesting 
result and it was assumed in a previous work [13]. 

Considering the above remarks and defining in the generating functional modification an 

alternative nonnegative constant C = 2<y2 J^ (which will be called condensate parameter) for 
further convenience, the expression ( 4.21| ) takes the form 



exp J I AA^^ 0*0 r> a (y) 



\4 , 



( 27 2 J 

For the above construction the vacuum in the zero momentum limit has the form 



= exp 



El A a+ A a+ -\- - A a+ A a+ -4- R a+ A L ' a+ -A- ir a+ r a+ 
9 0,1 ^0,1 ^ 9^0,2^0,2 + n ^0 + <c c 



a=l 



0). 



(4.22) 
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Then in the zero momentum limit this state has exactly the same form of the one proposed 



previously |14| , 15], the only difference appearing is introduced in the way in which the zero 



momentum limit is taken. That previous work [14, |l5j corresponds to the particular case in 
which a = 1, C3 (a = 1) = C\. 

The modification to the ghost generating functional of the free-particle Green functions 



is the same calculated previously [14, 15 1. It is because the ghost sector in the vacuum state 
remains unchanged. As was mentioned there, for the value C3 (a, 0) = 1 as was fixed here, there 
is no modification of the ghost sector and its generating functional and propagator remain the 
same of the usual PQCD. The consequence of this unchanged sector was mentioned in Section 
3.2, that is it guarantees the fulfillment of a Ward identity in the one loop calculation for the 
gluon self energy. 

Then, after all the previous analysis, it is concluded that the vacuum state introduced 
( 4.22; ) modifies the usual perturbation theory only through a change in the gluon generating 



functional of the free-particle Green functions in such a way that the total gluon propagator 
(including the usual perturbative piece) takes the form 



(2< 



^2 (fiV - (1 - a) ) - >('<) (k)(] lir 



ex.p{—ik (x — y)} . 

(4.23) 



in which the first term is the usual perturbative propagator. 

Finally it is possible to analyze the consequences of the present generalization on previous 
calculations, for (G 2 ), quark and gluon masses, done in Chapter 3. 

For the mean value of G 2 and the quarks masses, after performing calculations like the 
ones realized previously, was found that the results are exactly the same obtained before, 
independently of the a parameter values. However for the gluon mass is obtained the following 
result 



2 3(l + a)g 2 C 

m r = -. . 

(2tt) 4 

This result, a dependent, does not represent a contradiction because the gluon mass is 
tachyonic for all the admissible values of the a parameter (the non-negative ones), then it is 
not a physical magnitude and could depend on the gauge choice. However at the present state 
of the analysis the general proof of the gauge invariance of the theory is lacking and will be 
considered in future works. 
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Chapter 5 



Summary 



and has been calculated previously 
The condensate parameter g 2 C was determined to 



This work considers the extension of the former ones IS ]-[15] . The results obtained can be 
summarized as follows: 

1) The elements of the modified perturbative expansion for QCD initiated in were 
specified. The diagram technique embodies non-perturbative information associated to a non- 
trivial vacuum in which exists a condensate of soft gluon pairs, as found in [14, 15 1. 

2) The former calculation of the mean value of G 2 was completed by including all corrections 
of order g 2 , disregarding the terms not including the condensate effects. The result was the 
same obtained previously in the tree approximation 14 ] . The mean value of G 2 is a relevant 
quantity for the description of the QCD vacuum |17| 
considering non-perturbative methods, 
reproduce the accepted value of (g 2 G 2 ). 

3) The contribution of the condensate to the gluon self-energy, up to the g 2 order, was 
calculated. A tachyonic result for the gluon mass was obtained, for the parameter values 
employed in the construction of the initial state defined in [14, 15|. This because the condensate 
parameter C was determined as real and non-negative in that work. The polarization of the 
tachyonic mode is transverse. As remarked before there are recent arguments aiming at the 
possibility of the existence of a tachyonic gluon mass [26, 27 [, and its role in improving the 
quark antiquark potential in bound state models of mesons 27]. 

4) An evaluation for the quark masses, under the influence of the condensate, was done. 
Various solutions were obtained for the masses of these particles, among which only one is 
real for all the values the current mass parameter. Employing this particular solution some 
interesting predictions were derived. In particular values of near one third of the nucleon 
mass, for the light quark masses, were obtained. Further, the effective quark masses were 
evaluated as determined by the known current mass for each flavor. By mean of these values, 
the ground state energies of some hadron resonances, reported in |TlJ], were estimated through 
the simple addition of the masses of their known constituent quarks. It could be remarked 



34 



that the results obtained are reasonable good ones in comparison with the experimental data, 
after taking into account the high errors in the experimentally determined Lagrangian masses. 
This makes possible to sustain the hypothesis about the almost glue nature of the constituent 
mass values of the light quarks. The small contribution of the binding potential energy of such 
resonances is also indicated. This conclusion then supports the applicability of non-relativistic 
approximations in their study. 

5) The effective potential was determined, in the approximation considered in a previous 
work |l3|], as a function of (G 2 ). As it was obtained previously [jO]] the condensate is generated 
spontaneously from the zero field vacuum, as in the Savvidy Model. That is, it was shown that 
the effective action has a minimum for a mean value of G 2 different from zero. 

It is worth resuming also the main aspects of the modified expansion that have been clarified 
after the initial work [13|: 

A) The Feynman expansion, depending on the condensate parameter and a gauge parameter 
a = 1, corresponds to the Wick expansion around a physical state of the free QCD. 

B) Also in [14, [D| the condensation parameter C was defined as real and positive. This 
result determined the tachyonic value for the gluon mass evaluated here and the masses for up 
and down quarks as 1/3 of the nucleon mass. 

C) Another aspect which is important to underline is that the selection of the parameters 
in the initial state [14, 15 1 was designed also to impose the absence of a modification for 



the free ghost propagator. However, as it followed from the present analysis, this property 
in turns is related with the fact that a condensate modification of the ghost propagator can 
produce a longitudinal contribution to the self-energy. But, such a term should not exist for 
the transversality condition of the polarization operator (Ward identity) to be obeyed. Thus, 
its appearance could break manifestly the gauge invariance. Therefore, the present work is also 
given foundation to the non-modified ghost propagator choice considered in [13|-[15[. 

5) Finally, it was explored the generalization of the theory for arbitrary values of the a gauge 
parameter. A particular vacuum state for a Perturbative QCD was proposed, it is formed by 
a coherent superposition of zero momentum gluon pairs, and it is a generalization of the one 
proposed in a previous work . lt 

was analyzed that the proposed state satisfies the BRST 
physical conditions imposed by the operational quantization method for gauge fields, developed 



by Kugo and Ojima |2(J, 21]. The modifications in the perturbation theory introduced by the 



proposed vacuum state were calculated for arbitrary values of the a gauge parameter. The 
parameters introduced in the definition of the vacuum state were fixed in such a way that they 
reproduce the modification to the gluon propagator initially proposed and based in a previous 
work (l^] starting from a functional analysis. The condensation parameter C, that modifies 
the standard propagator multiplying a term containing a delta function in the four-momentum 
p, was determined to be real and nonnegative. The consequences of the generalized theory on 
previously calculated magnitudes were analyzed. For the calculation of the mean value of G 2 
and the quark masses was obtained that the results are identical to the ones determined for 
a = 1. The expression for the gluon mass was also calculated and it is determined to be gauge 
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dependent and tachyonic for all the values of the a parameter. The general proof of the gauge 
invariance of the theory is lacking and will be study later on. 

At the present stage some possible lines for the extension of the work can be visualized: 

1) To explore in more detail the various possibilities for the quantization of the zero modes 
of the gluon field. It could allow to bypass the employment of the zero momentum limit used 
in the construction of the initial state, and also perhaps permit to simplify its introduction and 
the study of alternatives for the modified expansion describing different physical ground states. 

2) Improve the study of the effective action as a functional of the condensate parameter, in 
order to search for a variant of the leading logarithm model useful for the investigation of field 
configurations associated to inter-quark strings, nucleon-nucleon potentials, etc. 

3) Start the introduction of the field and coupling constant renormalization in the scheme. 
This would allow to investigate if the infrared behavior of the running coupling of the PQCD 
can be regularized by the modified expansion. The evolution of the masses at high energy will 
be another task to be worth considering. At last the determination of the absolute value of the 
constant C would be also possible after finding the running coupling evolution. 

4) To investigate the possibility for a derivation of existing successful bounded state models 
for the heavy quark mesons [28]-|K]], as realized by the a ladder approximation for the Bethe- 
Salpeter equation within the proposed modified expansion. The presence of a tachyonic gluon 
propagator in the approach (which is argued to have the effect of introducing a linearly rising 
component in the inter-quark potential P7fl ) and the obtained constituent values for the light 
quark masses already support such a possibility. 
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Appendix A 

Diagrams for the evaluation of (G 2 ) 



For the exposition of diagrams appearing in the (G 2 ) calculation, only the gluon and quarks 
will be taken into account. The terms related with ghosts do not contribute. This conclusion 
can be extracted easily following the diagrams related with quarks, since the associated ghost 
diagrams are similar. All the contributions to (G 2 ) of diagrams not including the condensate 
parameter are disregarded, as was justified in the introduction. 

Expanding the exponential in the vertices in ( 2.10 ), up to order g 2 for the gluon and quark 
contributions it is obtained 



1 + ¥L^±± + S a b c d 6 5 5 5 



+- 



3H 2 5J a 5J b 5 J, 
5 



+ 



S? a j 5 5 



4!i3 5J a 5J b 5J c 5J d i 2 5^5J a 5{-Cj) 



5 Sdef 5 



5 5 



3\i 2 5J a 5J b 5J C 3\i 2 5J d 5J e 5J f 
5 5 S idj 5 5 



+ 



S?aj 5 5 



S?u 5 5 



+ 2 



3li 2 5J a 5J b 5J C i 2 5& 5J d 5 
where the norm (N) up to order g 2 is given by the expression 

N=l - 



i 2 S&SJaSi-b 



(A.l) 



^abcd $ 




5 5 5 



4!i3 5J a 5J b 5J c 5J d 



5 S def 5 



5 5 



3\i 2 5J a 5J b 5J C 3\i 2 5J d 5J e 5J f 



+ 



^iaj 5 



qQ c 



5ii5J a 5{-^)2\i 2 5i k 5J b 5{-^i) 



+2- 



Sg c 5 5 6 Sft 5 5 



id) 



i 2 5 J a 5 J b 5 J c 2H 2 S&5 J d 5 



Za[J,Z,Z] 



(A.2) 
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Now, the diagrams generated by the action of the functional differential operator in fl3.3| ) 
on (|A.1| ), after taking into account ( |A.2[) , are analysed. 

The third term in the expansion (|3.3[) is already of order g 2 . Up to this order also, it will 



only contribute with the first term in the expansion ( A.l ), which only generates the diagram 



shown in Fig. A.l 



F, x [CO] 



Figure A.l: Diagram 1. 



Substituting the expressions for propagators and vertices in the diagram of the Fig. A.l 
disregarding the terms in which the condensation parameter C do not appear and considering 
that the integrals of the form 



d D a 

I = I - — 1- = for a > (A.3) 

(-r) 



within dimensional regularization, the following result for the first diagram follows: 



576g 2 C 2 



d 4 x 



(A.4) 



where J d A x means the space-time volume. The symmetry factor of ( A.3 ) is F\ = — g. This 
diagram represents the tree approximation for the evaluation of the mean value of G 2 and it 



was calculated before in [13]. 



The second term in the expansion (3.3) contributes with the second and the fourth terms 
in the expansion ( A.l| ). 



With the second term of (A.l) two topologically independent diagrams appear which are 



represented in Fig A. 2 



F 2 X [^] + F 3 X [OO] 



Figure A. 2: Diagram 2. 
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With the fourth term of (A.l 



) , one diagram appears which is represented in Fig |A.3j . 



F 4 X [OO] 

Figure A. 3: Diagram 3. 



For the calculation of the first diagram in Fig. A. 2, when it is taken into account the ie 



prescription for the gluon propagator, is obtained a null contribution. 



The second diagram in Fig. A. 2 and the one in Fig. A. 2 generate a null contribution to 
the mean value of G 2 . This is due to the appearance of terms of the type f abc fi bc and T l ^' a S 1 ^ 
respectively which are equal to zero due to the antisymmetry of the structure constants f abc 
(similar arguments are used for the analogous diagrams having ghosts) and also due to the 
vanishing traces of the matrices T lha . 

Next, the first term in the expansion ( |3.3|) , up to the order g 2 , contributes with the first, 
third, fifth, sixth and seventh of the expansion ( |A.l ), as well as with the four terms of (A. 2). 

With the first term of (|A.1|), the diagram represented in Fig. A~~4 is obtained. 



R * [O] 



Figure A. 4: Diagram 4. 



With the third term of ( |A.1| ) and the first of ( |A.2| ), the diagram shown in Fig. A. 5 is 



obtained. It can be observed that the net effect of the normalization factor of the generating 
functional has been to cancel out the vacuum diagrams, the same occurs for the diagrams which 
follow. 

F 6 X [OO] 

Figure A. 5: Diagram 5. 

With the fifth term of ( |A.1| ) and the second of flA.2| ), the three diagrams represented in 
Fig. A. 6 appear; with the sixth term of ( A.l ) and the third of ( |A.2| ), the diagrams represented 
in Fig. |A.7| ar e obtained; finally for the action of (3.3) on the seventh term of (A.l) and the 
fourth of QA.2j) it follows Fig. 
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F 7 x [OO] 

+ F 8 X [Q~~~0] 



*[<Q] 



+ F 9 L . 

Figure A. 6: Diagram 6. 



F 10 x[O— 0]+F n x[<Q] 



Figure A. 7: Diagram for 7. 



F 12 X [OO] 

Figure A. 8: Diagram 8. 

Once, all the diagrams generated by the first term in the expansion ( |3.3| ) were examined, 
the results for their contributions are analyzed. 



For the diagram represented in Fig. [A.4| , in the momentum representation it is obtained: 

/ [ " ( 9tHtijk2 " S^g^ {-tC5 (k)) f d 4 x = 0, (A.5) 

Thanks to the antisymmetry property of the structure constants and the vanishing traces of 



the matrices T lJ,a , the first and second diagrams of Fig. A. 6, the first of Fig. A. 7 and the two 



of the Fig. A. 8 will also vanish. In this way, it is only needed to calculate the contributions of 



the following diagrams: the one in Fig. A.5, the third of Fig. A. 6 and the second of Fig. A. 7. 
After considering dimensional regularization and the ie prescription it was found a vanishing 
result for the contribution of that diagrams. 

It is important to remark that the same occurs for the diagrams generated by the ghost 
particles. Then the only diagram contributing to the mean value of G 2 , in the considered 



approximation, is the one in Fig. A.l. 
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Appendix B 

Longitudinal and Scalar Modes 
contribution 



The longitudinal and scalar modes contribution is determined by the expression 



(0 | exp \C* 3 (a, P) B a p A L A exp \ i I d 4 xJ^ a (x) (A^f + /* ^ (*) + 



+B a + [f; i}S}fl (x) + (1 - a) {A* (p ,x ) /* ^ (x) + B; fa x))] ) } 
x exp |i j d 4 xJ^ a (x) {Bi [f PiAlx {x) + (1 - a) (A fa x ) f n , Ljtl (x) + B^ fax))] + 

+A L ,ff Pi ^ (x)) } exp [C S («, P) B<i+ A^ + ] | 0). (B.l) 



Using the same methodology of the previous work ][L4|, 15], is first calculated 



exp U f d 4 xJ^ a (x) [f Pi ,s,» (x) + (1 - a) (A fa, x ) f Pi , L ^ (x) + B^ (r, x))] + 

+A^f Pi , L> , (x)) } exp [C 3 (a, P) BgJlg*] | 0) 
= exp |c 3 (a, P) (b£ - i J d 4 xJ^ a (x) f p . )LtlM (x)\ x (B.2) 

i [ d 4 xJ^ a (x) [f PiAlx (x) + (1 - a) (A fa,x ) f Pi , L ,„ (x) + B„ (pt,x))]) ) | 0), 



x ( A L J a+ - 1 



in the same way can be calculated the left part of ( |B,l[ ) 
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Substituting the right (|B.2|) and left parts in flB.l|), and introducing the following notation 



C* 
A+ 

a L 

h 
b-2 



= C* 3 (a,P), C = C 3 (a,P) 



A 



L,a+ 



A = Ab a , B + = B a +, B = B% 

Pi ' Pi ' Pi 



* J **J» (») [f^ (*) + (1 - a) (A* (p , , ) n ^ (*) + b; m)] , 
i I d A xJ^ a (x) [f Pi>s ,» (x) + (1 - a) (A (p , xq) f Pi ,L,„ (x) + B^ ( Pl ,x))] , 

i J d 4 xj^(x)f;^(x), 

i [ d i xJ^ a {x)f n ^{x), (B.3) 



the expression ( [B.lD takes the form 

(0 | exp \<0* (A + oi) (B + bx) } exp {c (s + + 6 2 ) (i+ + a 2 ) } | 0). (B.4) 

Following exactly the same procedure described previously [14, [Tq ], the following recurrence 
relation is obtained for (B.4). 



exp J C%6i + C (C*oi - o a ) (C*&i - 6 2 ) ^ \^C\ 2{2m) + |C| 



2(2m+l) 



(B.5) 



771=0 



(0 I exp {(Til? + C* n+1 C n+1 U<?bi -b 2 )A + (C*ai - a 2 ) B) } exp {c^+i+J | 0). 



To take the limit n — > oo it is considered |C| < 1, which implies the expression ( p.5[ ) takes 
the form 



exp < 



C*ai6i + C (C*oi - 02) (C*6i - 62) 



1 



1-|C| ! 



(0 I exp {(Til?} exp | 0). 



(B.6) 



Finally the notation (|B.3| ) is substituted in flB.q ), After that all the functions of f>i are 
expanded in the vicinity of pi = 0, keeping in mind that the sources are located in a space 
finite region, it is necessary to consider only the first terms in all the expansion. Then, for the 
expression ( |B.6| ), the result obtained is ( 4.1 6| ) . 
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